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We derive the thermomagnonic torque associated with smooth magnetic textures subjected to 
a temperature gradient, in the framework of the stochastic Landau-Lifshitz-Gilbert equation. Our 
approach captures on equal footing two distinct contributions: (1) A local entropic torque that 
is caused by a temperature dependence of the effective exchange field, the existence of which had 
been previously suggested based on numerics and (2) the well-known spin-transfer torque induced 
by thermally-induced magnon flow. The dissipative components of two torques have the same 
structure, following a common phenomenology, but opposite signs, with the twice larger entropic 
torque leading to a domain-wall motion toward the hotter region. We compare the efficiency of the 
torque-driven domain-wall motion with the recently proposed Brownian thermophoresis. 
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Introduction .—Spintronics aims at the active con¬ 
trol and manipulation of spin degrees of freedom in 
condensed-matter systems. 1 Spin dynamics can be in¬ 
duced by various means, e.g., an external magnetic fielcP 
or an electric current (in conducting systems, such as 
magnetic multilayers)! 3 A thermal flux offers another 
general instrument to excite magnetic dynamics, which 
forms the central subject of the field dubbed “spin 
caloritronics.”® The spin Seebeck effect, in which a tem¬ 
perature gradient generates spin current, is a prime ex¬ 
ample of a spin-caloritronic phenomenon, which has been 
observed in a ferromagnetic metalP a semiconductor, 6 ! 
and an insulator. 7 Thermal control of a magnetic system 
has an advantage over an electronic control that is nec¬ 
essarily associated with undesired energy dissipation due 
to electronic continuum. 

A ferromagnetic domain wall is a classic example of 
topologically-stable magnetic textures; it has been ex¬ 
tensively studied due to a fundamental interest as well 
as practical motivations exemplified by the racetrack 
memoryP Recent experiments have demonstrated that 
a temperature gradient drives domain walls to the hot¬ 
ter region! 9 On the theoretical side, a thermomagnonic 
torque on a domain wall or, more generally, a mag¬ 
netic texture has been addressed from several perspec¬ 
tives. Hinzke and Nowak® predicted a thermally- 
induced motion of a domain wall by numerical simu¬ 
lations based on the stochastic Landau-Lifshitz-Gilbert 
(LLG) equation; 11 Schlickeiser et a/.® argued later that 
the motion is driven by the so-called entropic torque that 
originates from the temperature dependence of the ex¬ 
change stiffness. Yan et al . 13 derived an adiabatic spin- 
transfer torque (STT) on the domain wall exerted by 
magnons based on the conservation of angular momen¬ 
tum: Thermal magnons adjust their spins toward the lo¬ 
cal magnetic order and react by exerting a torque on the 
domain wall. Kovalev and Tserkovnyak SI constructed a 
phenomenological theory for the thermomagnonic torque 
that includes the so-called “/3-type” dissipative correction 
to the STT, which has been further elaborated analyti¬ 
cally in the framework of the stochastic LLG equation.® 


Domain Wall Motion 



FIG. 1. (Color online) An illustration of the motion of a do¬ 
main wall in a one-dimensional ferromagnet under thermal 
bias. The unit vector n (blue arrow) in the direction of the 
local spin density precesses around the equilibrium position 
due to thermal magnons. As the temperature increases, the 
length of the time-averaged n (yellow arrow) decreases, and 
thus the effective exchange field decreases as well. The tem¬ 
perature dependence of the effective field engenders a torque 
on the domain wall, pushing it to the hotter region. The STT 
acting in the opposite direction is not strong enough to reverse 
the direction of the motion within the LLG treatment. 


In this work, we reconcile these theories (resolving 
some of the conflicting findings) for the thermomagnonic 
torque within the LLG phenomenology (applicable at 
temperatures T « T c , the Curie temperature), which 
serves as their common underlying framework. In the 
following, we start by formulating the model and sum¬ 
marizing our main results: The first is the analytical 
identification of a contribution to the thermomagnonic 
torque that is caused by temperature dependence of the 
effective exchange field, the existence of which had been 
previously supported by numerics.® The second is the 
expression of the total thermomagnonic torque in terms 
of the magnon current density and the gradient of the 
magnon-mimber density, which point to the physical ori¬ 
gin of each term. We then provide the detailed deriva¬ 
tions of the results. The thermomagnonic torque adds 
an additional term in the equations of motion for the 
collective coordinates, which parametrize slow modes of 
a magnetic texture. We study the motion of a domain 
wall (sketched in Fig. [TJ as an example. Lastly we dis¬ 
cuss other possible effects on thermal motion of magnetic 
domain walls. 
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Dynamics of a ferromagnet at finite temperature can 
be described by the stochastic LLG equation: 11 


r st = — diJf. The spin current can be split into the lon¬ 
gitudinal and transverse (relative to n^ 0 ^) components: 


5(1 + anx)n = n x (h + h th ), (1) 


where n is the unit vector in the direction of a local spin 
angular-momentum density s = sn, a is the Gilbert¬ 
damping constant^ 6 h = —dU/dn is the effective field 
conjugate to n. For the ferromagnet’s energy, we take, 
for concreteness, U[ n] = J dV[A\Vn\ 2 + KijUiUj]/2 , 
where temperature-independent constants A > 0 and 
Kij = Kji parametrize the exchange stiffness and mag¬ 
netic anisotropies, respectively. Fluctuations of spins at 
finite temperature are modeled by the Langevin field 
h th ; the fluctuation-dissipation theorem 17 relates the 
Langevin field to the damping term, which manifests in 
the correlator of h th (setting ks = 1): 


(h lh (r, cj)h lh (r / , a/)) 


27 rSijashuj 

tanh(7kj/2T) 


r')5(u; — uj'). 

( 2 ) 


Main results .—We split n entering Eq. 0 into two 
orthogonal components: n = a/1 — £n 2 n(°) + 5n. The 
slowly-varying component n^ 0 ) describes a smooth mag¬ 
netic texture; the fast component (in represents a small 
deviation of the local spin density from that is caused 
by thermal agitation. Thermal magnons constitute the 
fast component 5n, whose effect on the magnetic texture 
is captured by the thermomagnonic torque in the (tem¬ 
porally) coarse-grained LLG equation. 

Averaging over the quadratic terms in the fast compo¬ 
nent 5 n yields the low-temperature LLG equation (pro¬ 
jected transverse to n^ 0 ^) for the slow component n^ 0 ^: 


s(l + 1 


' A JLL' 


- An (0) x diU^di(Sn 2 ) 
+ A(Sn x V 2 (in), 


in the exchange-magnon approximation, within which 
only the exchange energy is retained for thermal 
magnons. Here, h^ 0 ) = iV 2 n^ — is the effec¬ 

tive field for at zero temperature.^ 


The second term (3b) on the right-hand side is a 


magnonic torque that constitutes our first key result: 


= -(2ftA/s)n<°> x (dipdi) 


,(o) 


( 4 ) 


which exists due to temperature dependence of the ef¬ 
fective exchange field AV 2 (V 1 — £n 2 n(°)) through the 
magnon-number density p = s(Sn 2 )/2h. This contribu¬ 
tion is analogous to the so-called entropic torque, which 
has been argued in Ref. [12] to govern the domain-wall 
motion. 


The third term (3c), which is identified as the 
magnonic STT, 13 S can be expressed in terms of the 
divergence of the spin current Jf = — A($n x did n): 


J? A[n^ • (5n x di5n)]nW 

-V°) x (5n(^n(°) .<Sn)). l >) 

This finally leads to the following expression for the 
magnonic STT: 

r st m h (J • V)n (0) + (hA/s) n (0) x (dipdi) n (0) , ( 6 ) 

to first order in the spatial derivative of the texture. 
Here, Ji m (A/h)[n^ • (dn x did n)] is the magnon-flux 
density evaluated in the absence of a (slow) magnetic tex¬ 
ture. For the circular exchange magnons, we similarly 
evaluated (dn(3in( 0 ) • Sn)) = (hp/s)diU^°\ 

Adding the two torque contributions, Eqs. 0 and ([ 6 ]), 
yields the total thermomagnonic torque r = r ex + r. 

t = H (J • V)n(°) - (SA/«)n(°) x (d iP di) n<°), (7) 

which is a central result of the paper. The benefit of 
using macroscopic hydrodynamic variables J and Vp for 
the thermomagnonic torque is twofold. First, their di¬ 
rections relative to a temperature gradient are apparent: 
J = —cjVT and Vp = c p VT, with positive constants 
cj and Cp. As a result, the relative sign of two terms 
in Eq. ([7| for r is unambiguously established, which was 
not possible in the phenomenological approach.^Second, 
it sheds light on the physical origin of the thermally- 
induced magnetic-texture dynamics. In particular, a fer¬ 
romagnetic domain wall is driven to the hotter region at 
the speed v = HAdip/as 2 (normal to the domain-wall 
orientation), which is governed by the gradient of the 
magnon-number density. 

The solution for the fast component of the stochastic 
LLG equation, which we obtain below by the standard 
approach} 15 ® yields in d dimensions: 

_ h 1 _ did s 

CJ ~ 2 d ir d / 2 ahX d ~ 2 ’ Cp ~ 2 d + 1 it d / 2 AhX d ~ 2 ’ [ ’ 

where A = xJhA/sT is the thermal-magnon wavelength. 
Here, I d = [1/T(1 + d/2)} f 0 °° dp r] d / 2 e^{e^ - l)~ 2 is a 
numerical constant (which would need to be regularized 
at low energies for d < 2). For d > 2, it is given by 
the Riemann zeta function,^ I d = C(d/2), with 1 3 ~ 
2.612 for the physically most relevant case. Note that the 
first term oc J in the thermomagnonic torque r, Eq. 0, 
becomes dominant in the limit of damping a -» 0 as then 
cj oc while c p remains finite. 

The thermomagnonic torque ([T]) is composed of two 
terms linear in magnetic texture: The one associated 
with the magnon-flux density J is the reactive component 
of the adiabatic torque ! 13 14 The other, caused by the gra¬ 
dient of the number density p, is known as the u /3- typ e” 
dissipative contribution to the adiabatic torque ! 14 15 It 
has been introduced in the form of /3hn^ x (J ■ V)n(°) 
by the phenomenological approach . 14 Our derivation in 
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the framework of the stochastic LLG equation yields 
/3J = — (A/s)Vp, where /? = [A/ s)c p /cj = (d/2)<a.®l 
Note that the overall sign of /? is positive while the STT 
contribution © alone would result in a negative / 3.®1 
Thermal magnons. —In order to calculate the magnon 
flux, J, and number, p, densities entering in Eq. 0 , we 
have to solve the stochastic LLG equation 0 for the fast 
component: 

s(5n + x Sn) = n^ 0 ) x AX 2 Sn + n^ 0 ) x h th , (9) 

to the linear order in h th and supposing uniform constant 
order n^ 0 ). When the average thermal energy of magnons 
is much larger than the magnon gap, T hK/s , which is 
usually determined by the strength of the anisotropy if, 
it suffices to retain only the exchange energy of magnons. 
The fast mode $n can be combined into one complex field 
= 5n • (ei + ie 2 ) within the local orthonormal frame 
{ei,e 2 ,e 3 = n^ 0 )}. The complex field d satisfies the 
stochastic Schrodinger equation with dissipation : 15 

is[ 1 + ia)ip = —AX 2 ^ — h. (10) 

Here the complex Langevin field h = h th • (ei + ^ 2 ) has 
the correlator, (h*(r, t)h(r) t')} = 2 (h th (r, t) • h th (r',t')) 
(retaining only the transverse components h th _L n^ 0 )). 
The magnon-flux density, J, and the gradient of the 
magnon-number density, Vp, are given by 

J = (A/fi)Im(^*V^), Vp = (s/fi)Re(^*V^), (11) 

which can be easily obtained from Eq. (Jl0|. 
Fourier transforming along the transverse coordi¬ 
nates r_L = (y,z) and the time axis, ^(x, rj_, t) = 
f [d d ~ 1 k±duj/(27r) d ]e zk - L ' rj -~ 1 ' ujt 'i/j(x, kj_, lj), yields the 
Helmholtz equation: 

A(d 2 T k 2 )d(x, k_L, u) = —h(x, kj_, c«j), (12) 

with the complex-valued wave vector k 2 = scj (1 + 
id) IA — fcj_. The stochastic field correlator is given by 
(h*[x, kj_, oo)h[x) k(j_, w r ) oc 2[27r) d S(k± — Id, ), with the 
remaining factors carrying over from Eq. pi. Employ¬ 
ing the one-dimensional Green’s function G[x — y) = 
iex.p[ik\x — y\)/2k [with Im(fc) > 0], we finally obtain 
after straightforward algebra 

(0‘(r,t)ViKr.t» = (f -i) 2 4t ,/ 2 < 13 > 


Dynamics of magnetic textures .—The deterministic 
magnon-averaged LLG equation for a general smooth 
texture is given by 

s(l + an(°) x )h(°) = n(°) x + r, (14) 


in terms of the full thermomagnonic torque r = r ex +r st , 
Eq. 0 . The theory can be further simplified by focus¬ 
ing on the relevant collective coordinates parametrizing 
dynamics of interest. A canonical example of this is 
a quasi-one-dimensional domain-wall motion described 
by the azimuthal angle <L(£) at the wall’s center with 
position X(t)P In general, encoding the dynamics of 
a texture id 0 ) (r, t) with a set of coordinates, q (t) = 
{qi(t), q2{t), ••• }, the motion of rd°) reflects evolution 
of the coordinates: rd 0 ) = (pdrd°)/ckp. Taking the inner 
product of Eq. (14) with id 0 ) x (<9rd°)/c%) followed by 
spatial integration yields the equations of motion for the 
collective coordinates: 


Gq-Tq + F + F m = 0, (15) 


showing the interplay between gyrotropic, viscous, con¬ 
servative, and magnon-induced forces acting on a slowly 
varying texture.®! Here, Gij = s f dV rd°) • (c^.rd 0 ) x 
c^rd 0 )) is the antisymmetric gyrotropic tensor®! T^ = 
as f dV(d qi d°) - d Qj d 0 )) is the symmetric viscous tensor; 
Fi = —d qi U is the conservative force; and, finally, 


FA 1 = 


h J dV[ n (0) x (dn (0) /%)] ' (J ' V)n (0) 

(. HA/s ) J dV(dn^/dqi) ■ (djpdj)^ 0 


(16) 


is the thermal-magnon induced force. 

As an illustrative example, we now specialize to 
the case of a domain-wall motion under the temper¬ 
ature gradient VT = c^Tx, which is illustrated in 
Fig- HI For the easy-:r;;-plarie easy- 2 -axis ferromagnet, 
with energy U[ n] = f dV(A\Vn\ 2 + K y n 2 y - K z n 2 z )/2, 
K y ,K z > 0 , a domain wall is a topologically-stable 
equilibrium defect. It is explicitly given by rd°) = 
(sin# cos <L, sin# sin <L, cos#) with 4> = 0 and cos# = 
tanh[(x — X)/A], where A = yJ~AjKf z is the width of 
the wall. Slow dynamics of the domain wall are well de¬ 
scribed by two collective coordinates: the position X and 
the azimuthal angle 4>.^The equations of motion (15) for 
X and <f> are 


keeping terms that remain in the limit a —>> 0 . Sub¬ 
stituting this result into Eqs. (11), we directly read out 
coefficients p). 

As an additional check, Vp can be alternatively 
calculated by considering magnons in the local ther¬ 
mal equilibrium at temperature T(r). Using the 
Bose-Einstein distribution function, the number den¬ 
sity of the thermal exchange magnons is given by p = 
f[d d k/(27rY)(e hAk / sT —1) _1 , whose gradient reproduces 
c p in Eq. p|. 


s[X - aA4>) = -hJ x + [Ky/ 2) sin(24>), (17) 

s(A<i> + aX) = ( hA/s)d x p . (18) 


The steady-state velocity of the wall below the Walker 
breakdown 2 ’ (in the three-dimensional case when wire’s 
lateral dimensions are larger than the thermal-magnon 
wavelength) is 


HA = 3C(3/2) d x T 
as 2 x ^ 167 T 3 / 2 asA 


( 19 ) 
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The associated steady-state azimuthal angle is given (at 
small biases) by 4> « ((3/2)d x T/167r 3 / 1 2 * 4 * * 7 a\K y . For a nu¬ 
merical estimate, taking s = h/ nm 3 , a = 10 -2 , A = 10 
nm, and d x T = 10 K/mm, the resultant velocity of 
the wall is v ~ 1.15 mm/s. The domain wall moves 
to the hotter region, with the driving force proportional 
to the gradient of the magnon-mimber density, which is 
consistent!^ with recent numerical findings. 12 

Discussion .—In developing the long-wavelength theory 
for magnetic-texture dynamics subject to thermal gra¬ 
dients, we have tacitly assumed that the thermal wave¬ 
length of magnons is much shorter than the characteristic 
length scale of the slow texture. This allows us to focus 
on the adiabatic torques that are first order in spatial- 
texture gradients. Further simplifications are afforded by 
assuming energy-scale hierarchy hK/s «T «T C , such 
that thermal magnons are dilute and dominated by ex¬ 
change interactions (hence circular). A strong anisotropy 
on the atomistic scale, however, would require one to ac¬ 
count for the noncircular character of magnons beyond 
the exchange approximation. Topological defects, such 
as domain walls, in this case, would, furthermore, be as¬ 
sociated with textures that can be sharp compared to the 
lattice constant. This would lead to strong nonadiabatic 
effects in magnon transport and torqu es, e ngendered by 
magnon reflection on sharp textures,^ 27 ^ that are be¬ 
yond our formalism. High temperatures approaching T c , 
furthermore, would invalidate the Landau-Lifshitz phe¬ 
nomenology as a starting point, calling for alternative 
approaches. 


Whereas we focused on deterministic texture dynamics 
in the long-wavelength limit of the bulk, a small topolog¬ 
ical soliton, e.g., a domain wall in a narrow wire, should 
generally behave as a particle immersed in a viscous 
medium. As such, it must exhibit Brownian motion due 
to random forces required according to the fluctuation- 
dissipation theorem. We have recently pointed out that 
a small antiferromagnetic soliton could drift to a colder 
region by the Brownian motion under a temperature 
gradient.Naturally, small ferromagnetic solitons be¬ 
have similarly. For a ferromagnetic domain wall of the 
type considered above, the drift velocity by the Brown¬ 
ian motion is vb = —(A/2 asa)d x T, where a is the mag¬ 
netic wire’s cross section in the yz plane. The ratio be¬ 
tween the deterministic torque-induced domain-wall ve¬ 
locity v, Eq. ( p~9] ) , and the stochastic Brownian drift vb is 
v/v B ~ — cr/AA. The thermomagnonic torque thus dom¬ 
inates for larger cross sections a (supposing rigid motion) 
and/or at higher temperatures, corresponding to shorter 
thermal-magnon wavelengths A oc T -1 / 2 . 
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